Ab-initio calculation method for charged slab systems 
using field-induced gaussian sheet 
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A new repeated-slab calculation method is developed to simulate the electronic structures of 
charged surfaces by arranging density-variable charged sheets in vacuum regions to realize a constant 
potential on the charged sheets and maintain the charge neutrality condition. The charged sheets are 
fabricated so as to screen an electric field from charged slabs; consequently, they act like a counter 
electrode composed of flat perfect conductors, modeling a tip of a scanning tunneling microscope 
or a reference electrode in an electrochemical cell. This method has the advantages of ease of 
implementation into existing repeated-slab programs and low computational costs. The availability 
of the method is demonstrated by applying it to a charged H2 molecule and charged Al(lll) metal 
and Si( 11 1) semiconductor surfaces. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 

A surface becomes an interesting stage for material 
transformations such as chemical reactions when the sur- 
face is placed under an external electric field and has pos- 
itive/negative charges. In the scanning tunneling micro- 
scope (STM), for example, the applied sample bias pro- 
duces an electric field between the surface and the STM 
tipi Such electric field often induces a variety of surface 
reconstructions, 2 desorption^ and the diffusion/rotation 
of adsorbatesi^ On the other hand, in electrochemistry, 
the manipulation of electric fields is a key factor in con- 
trolling chemical reactions*^ 

Ab-initio calculation using a repeated slab has greatly 
contributed to the exploration of the geometric and elec- 
tronic structures of a number of surfaces. Such calcula- 
tion employs a periodic boundary condition (PBC) per- 
pendicular to the slab; thus, it receives the benefit of nu- 
merous computational techniques developed in the calcu- 
lations of periodic bulk systems. However, most of these 
calculations are concerned with neutral surfaces, because 
various physical quantities such as energy and potential 
diverge to infinity unless charge neutrality is maintained 
within a unit cell in an electronic structure calculation 
based on the PBC schemed 

To treat charge-neutral surfaces in an external elec- 
tric field, Neugebauer and Scheffler 9 introduced artificial 
dipole sheets in vacuum regions of a repeated unit cell. 
Arranging the positive and negative charged sheets on 
the front and back sides of the slab, respectively, one 
can easily impose a constant electric field on the slab 
system. The calculation method for a charged slab was 
performed by Fu and Hoii^ In their method, artificially 
charged sheets, which have the same amount but an op- 
posite sign of charge to that of the slab, are placed on 
both sides of the slab. These artificially charged sheets 



terminate the electric field originating from the charged 
slab. All these computational recipes dexterously realize 
the charge neutrality and compensate the electric field in 
a unit cell, and they have been utilized in several studies 
for the elucidation of field effects on surfaces pliiSii&ii 

However, since these methods use sheets with fixed 
charge distributions, typically of plane Gaussian form, 
the electric field from a charged slab is in general not per- 
pendicular to the sheet and the potential shows a mod- 
ulation along the sheet. This feature becomes marked 
and produces some inconvenience when the surface of 
the slab is bumpy and the excess-charge distribution on 
the slab is far from uniform. For example, unless the 
thickness of the vacuum region is sufficiently large, the 
potential modulation on the sheet often produces unnec- 
essary electrostatic interactions between repeated slabs, 
which makes it complicated to analyze the energetics on 
the slab. Furthermore, in the case of a STM experi- 
ment or an electrochemical cell, since a STM tip and a 
reference electrode are normally made of conducting ma- 
terials, they have a constant potential but a nonuniform 
charge distribution, rather than a fixed uniform charge 
distribution (see Fig. 1(a)). To simulate the situations 
in which such a counter conducting material approaches 
a surface, it is desirable to adopt another calculational 
recipe using constant-potential charged sheets. 

Recently, Otani and Suginoi^ have developed a gen- 
eral method of simulating charged surfaces surrounded 
by screening mediums. They ingeniously combined a 
repeated-slab calculation with a real-space potential cal- 
culation to respond to any non-PBCs by releasing the 
charge-neutrality constraint of a unit cell. However, their 
method needs the regeneration of a potential in real space 
and requires additional computational costs. 

The purpose of this work is to demonstrate a new 
computational recipe for simulating a charged surface lo- 
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cated near a counter metallic electrode maintaining the 
charge-neutrality constraint by only inserting a constant- 
potential charged sheet within the repeated-slab scheme. 
Our approach provides the advantages of ease of imple- 
mentation into existing ab-initio calculation programs 
based on PBC and low additional computational costs 
for the charged slab calculation. The rest of this pa- 
per is organized as follows: in iflTl we explain how to 
construct a constant-potential charged sheet and then 
present formulas to realize it using an existing ab-initio 
calculation program. In AIIII we demonstrate how the 
present method works, by showing results for a charged 
H 2 molecule, charged Al(lll) and Si(lll) surfaces near 
an electrode. An attention with respect to the width of 
charged sheets is also presented. We give a brief sum- 
mary of this paper in flIVI 



II. FORMALISM 

In this section, we consider how to construct a charged 
sheet that has a constant potential on the sheet. When 
one connects a conducting material and a surface with 
an electric battery, counter charge is induced on the con- 
ducting material by an electric held from the charged 
surface, and this realizes a constant potential on a con- 
ducting material, as shown in Fig. l(a)ii& To construct 
the induced counter charge, we hrst obtain the electric 
field from a charged slab. 

We prepare an open space that is periodic in both x 
and y directions but nonperiodic in z direction, where 
the x, y and z directions correspond to the surface lateral 
and normal directions, respectively. The range of the z 
direction is limited to < z < 2L, and the center of the 
charged slab is located at z = L, as shown in Fig. ^b). 
We assume that the slab has reflection symmetry at z = 
L, for simplicity. The charged slab produces an electric 
field outside the slab, i.e., in the vacuum regions. The 
z component of such an electric field at the boundaries, 
z* = or z* = 2L, is obtained as 

E z (z*,g) = - d z G(z = z*,z',g)p slb (z',g)dz\ (1) 
Jo 

where we performed a Fourier transformation along the 
x — y direction and used the reciprocal component, g. 
Paib{z, g) is the electronic density of electrons and ions in 
the slab region. G(z, z',g) is the propagating Green's 
function obtained by solving the Fourier-transformed 
Poisson equation, (d z — g 2 )G(z, z',g) = —4:TtS(z — z'), 
where g = |g|. If we set the electronic potential to zero 
at the cell boundaries, G(z = z* , z', g) = 0, this Green's 
function can be expressed as 

G(z,z>, g ) = ^e-s\*-*'\ + ^ 
9 9 

e- 2Lg cosh(g(z - z')) - cosh(g(z + z' - 2L)) 



sinh(2L<?) 



Then, we arrange the charged sheets at the boundaries, 
z* = and z* = 2L. If the charge density on the sheet 
is determined so as to compensate E z {z*,g) 1 a constant 
potential is realized at the cell boundaries and the charge 
neutrality is naturally satisfied in < z < 2L on the basis 
of Gauss's theorem. Therefore, we employ E z (z*,g) to 
determine the charge densities on the sheets and adopt 
the half-Gaussian forms for the charge distribution along 
the z direction. 

Since the above-mentioned cell, which is composed of 
the charged slab and counter charged sheets at the cell 
boundaries, is neutral and produces no electric field out- 
side the cell, we stack such cells periodically along the 
z direction. In this case, the repeated counter charge 
density, p F (z 1 g), is given by 
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where 9 and a are the step function and the width of a 
Gaussian charge, respectively. A is a normalization con- 
stant of a half-Gaussian distribution, A = J °° e~^^ = 

^a. Inserting eqs. and into eq. © and Fourier- 
transforming the result along the z direction, we obtain 
for g ^ 
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where k is the reciprocal component along the z direc- 
tion and (f> s ib(k,g) is the Coulomb potential due to the 
slab charge, 47r ^^gf ■ no represents the excess electron 
number in a slab given by no = ^ unit p s ib{v)dY. Posi- 
tive and negative no's correspond to an excess and a de- 
ficiency of electrons in the slab, respectively. Hereafter, 
we call this type of counter charged sheet a "field-induced 
Gaussian-charge (FIGC) sheet". 

In the repeated-slab calculation, the total charge den- 
sity is obtained as the sum of the charge density of the 
FIGC sheet, /o_p(fc,g), and that of the slab, p s ib{k,g). 
Such a sum enables the electronic-structure calculation of 
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FIG. 1: (a) Schematic picture of charged surface system considered in this work. The surface and the counter electrode are 
connected with a cathode and an anode of an electric battery, respectively. Arrows denote an electronic field, E, in the vacuum 
region, (b) A slab setup employed in the present calculation to simulate the charged surface shown in (a). The center of the 
slab is located at z = L and each surface of the slab has a charge +a. Two counter charged sheets are set at both boundaries 
and each electrode has the counter charge — a. 



charged surfaces located nearest to a constant-potential 
electrode by a standard ab-initio method based on the 
PBC. By only adding pp(k,g), we can obtain physical 
quantities such as the electron distribution, the electronic 
potential, forces on ions, and the total energy of the sys- 
tem composed of the charged surface with the contour 
electrode^ 



III. RESULTS AND DISCUSSION 

In this section, we demonstrate how the present scheme 
works. We consider a charged H2 molecule, and charged 
Al(lll) and Si(lll) 1 x 1 surfaces. These systems are 
arranged at the center of a repeated unit cell and the 
atomic positions of all systems are optimized. We per- 
form a standard ab-initio calculation based on the den- 
sity functional theory (DFT) 18 - 19 using a plane-wave ba- 
sis set with an energy less than 25 Ryd. for the H2 sys- 
tem and 16 Ryd. for the Al(lll) and Si(lll) systems. 
The present scheme is implemented into the computa- 
tional code "Tokyo ab-initio program package" (Tapp). 20 
The Perdew-Burke-Ernzerhof (PBE) 21 functional in a 
generalized gradient approximation (GGA) is employed 
as an exchange-correlation potential for the H2 system, 
while the Perdew-Wang 2 - functional in a local density 
approximation (LDA) is adopted for the Al(lll) and 
Si(lll) systems. The ionic cores are described by ul- 
trasoft pseudopotentialsS in all cases. 

Figure |2fa) shows the setup of a H 2 molecule in a 
unit cell. Unit cells having 10.0 x 10.0 x 12.5 and 
10.0 x 10.0 x 18.75 a.u. 3 sizes are adopted. The H 2 
molecule is surrounded by vacuum in all directions. The 



counter charged sheets are set at unit cell boundaries of 
the z direction. We ionize this molecule positively as 
H^ 01 by extracting 0.1 electrons. Figures^b) and[3{c) 
show the counter maps of induced charge density cal- 
culated at the cell boundaries, when the charged sheets 
are placed 5.5 and 8.7 a.u. apart from the nearest hy- 
drogen atom of the molecule, respectively. In Fig. Hfc), 
since the distance between the molecule and the counter 
charged sheet is large, the electric field generated by the 
molecule uniformly touches the sheet and the charge dis- 
tribution on the sheet becomes uniform. However, when 
the sheet approaches the molecule, reflecting the nonuni- 
form charge distribution of Hj , the induced charge 
shows a spacial modulation, as shown in Fig. Hfb). 

Next, we consider positively charged Al(lll) and 
Si(lll) surfaces. FiguresEJa) and^b) show the excess- 
charge distributions of these systems viewed from the 
[101] and [110] directions, respectively. Both surfaces are 
located about 10 a.u. from the counter electrodes. The 
excess-charge distribution is calculated as the difference 
in the charge density between charged and neutral sys- 
tems. In the case of the Al(lll) metal surface, since the 
electrons are weakly bound to atoms in metals, the excess 
charge is uniformly distributed along the surface. As a 
result, the counter charge on the sheet becomes uniform 
an parallel to the surface. On the other hand, in the case 
of the Si(lll) surface, because dangling bonds have the 
highest energy in the surface, the inserted positive charge 
is mainly localized around dangling bonds of surface Si 
atoms. Thus, the counter charge on the sheet focuses at 
positions just above and below the dangling bonds. 

Then, we consider how a constant zero potential is 
realized on the counter charged sheet. In the present 
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FIG. 2: (a); Schematic view of charged H^ 0,1 molecule arranged in repeated unit cell. FIGC indicates the position of field- 
induced Gaussian-charge sheets, which are located at both the top and bottom of the cell, (b) and (c) are contour maps of 
charge distribution on the sheets when the sheets are located (b) 5.5 and (c) 8.7 a.u. from a nearest hydrogen atom. Denser 
color corresponds to higher density. 
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FIG. 3: Cross-sectional [101] and [110] views of excess-charge 
distributions for (a) Al +01 (111) and (b) Si +01 (111) surface 
slabs. The charge distributions of the counter charged sheets 
are also shown. Blue and red areas correspond to positive and 
negative charges, respectively. Solid circles indicate atoms on 
the display plane, while open circles indicate atoms out of the 
plane. A Gaussian-charge width of a — 1.0 a.u. is adopted in 
both (a) and (b). 



repeated-slab scheme, there exists one artificial param- 
eter: the width of the Gaussian-charge sheet, denoted 
as a in cq. J3J). This width is introduced to decrease 
the number of the plane-wave basis set and, thus, the 
computational time. However, owing to a finite a, the 
zero-potential condition assumed in the present method 



is slightly released. This situation can be understood 
by observing the potentials on the sheet. The poten- 
tial at z = is, for example, obtained by the sum of 
the fc-components of the total potential as cf>(g, z = 0) = 
X]/bL-oo <HS' k)/2L. Using a simple calculation, we ob- 
tain 

0(g,« = O) = — <J>slb(g,k') 

k' — — OQ 

OO ( ka \2 

x{l-|tanh(L.g) £ (6) 

At a = 0, it can be analytically derived that the 
right hand of eq.© equals zero using the series for- 

mula E^-oo x^Tn? = xtJLfrx) ' where X is an ar " 
bitrary real number and n is an integer. Hence, the zero- 
potential condition is satisfied on the sheet. 

On the other hand, in the case of a finite a, we show in 
Fig. 0] the calculated (f)(g,z = 0) for the charged H^ 01 
system shown in Fig. I^a). Here, the charged sheets are 
arranged 5.5 a.u. apart from the hydrogen molecule. The 
4>{s, z = 0) values using the conventional fixed Gaussian- 
charge-sheet method are also shown for reference. This 
0(g, z = 0) represents the spacial modulation of potential 
on the counter charged sheets, corresponding to the wave 
vector g. The result using the fixed Gaussian-charge 
sheet shows a potential modulation of about 20 meV, 
thus, not realizing the constant zero potential. On the 
other hand, in the case of the present scheme, it is seen 
that the potential modulation of each wave vector is sup- 
pressed and the potential modulation decreases with a. 
Therefore, a zero potential can be realized as long as a 
small a is employed. 



IV. SUMMARY 

To calculate the electronic structures of charged sur- 
faces located near a flat and metallic counter electrode, 
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we developed a new method that uses neutralizing vari- 
able charged sheets to realize a constant potential on the 
electrode. Since the present method is constructed on 
the basis of the PBC and a neutral unit cell, this method 
has the significant advantages that the implementation 
into existing repeated-slab programs is easy and the ad- 
ditional computational time is short. We have demon- 
strated that the present scheme works well not only for 
metal surfaces but also for semiconductor and insulator 
surfaces, and surfaces with adsorbates. We hope that the 
present scheme will be implemented in various programs 
and used for the study of charged surfaces. 



FIG. 4: Fourier components of total potential on counter 
charged sheet at z — 0, 4>(g, z — 0) , for various in-plane wave 
vectors, g. A charged H} molecule is located at the center 
of a unit cell, 5.5 a.u. apart from the sheets. We adopt g char- 
acterized by the an integer m, as g = go (m, m), where go is 
the length of the in-plane fundamental reciprocal vector. Solid 
and open circles indicate the potential components obtained 
by the FIGC scheme using the widths a = 0.3 and a = 0.5 
a.u., respectively, whereas open squares denote those obtained 
by the conventional fixed Gaussian-charge-sheet method us- 
ing a width of a — 0.5 a.u.. 
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